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ABSTRACT 

The W3 GMC is a prime target for investigating the formation of high-mass stars and clusters. This second 
study of W3 within the HOBYS Key Program provides a comparative analysis of subfields within W3 to further 
constrain the processes leading to the observed structures and stellar population. Probability density functions 
(PDFs) and cumulative mass distributions (CMDs) were created from dust column density maps, quantified as 
extinction Ay. The shape of the PDF, typically represented with a lognormal function at low Ay “breaking” to a 
power-law tail at high Ay, is influenced by various processes including turbulence and self-gravity. The breaks 
can also be identified, often more readily, in the CMDs. The PDF break from lognormal (Ay(SF) ~ 6 - 10 mag) 
appears to shift to higher Ay by stellar feedback, so that high-mass star-forming regions tend to have higher 
PDF breaks. A second break at Ay > 50 mag traces structures formed or influenced by a dynamic process. 
Because such a process has been suggested to drive high-mass star formation in W3, this second break might 
then identify regions with potential for hosting high-mass stars/clusters. Stellar feedback appears to be a major 
mechanism driving the local evolution and state of regions within W3. A high initial star formation efficiency 
in a dense medium could result in a self-enhancing process, leading to more compression and favorable star- 
formation conditions (e.g., colliding flows), a richer stellar content, and massive stars. This scenario would be 
compatible with the “convergent constructive feedback” model introduced in our previous Herschel study. 

Subject headings: ISM: dust, extinction - ISM: individual (Westerhout 3) - Infrared: stars - Stars: formation - 
Stars: early-type 


1. INTRODUCTION 

The Giant MolecularCloud (GMC) W3 is rich in high-mass 
star activity (e.g. jMegeath et al.l2008l) and its relatively prox- 
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imity (~ 2kpc; e.g. , lHachisuka et al.1 120041 IXu et all 120061 
iNavarete et al.ll201 1) makes it a prime target for the study 
of cluster and high-mass star formation. W3 co ntain s high- 
mass stars in various evolutionary stages (e.g. jTieftrunk et all 
119971) . The eastern high density layer (HDL) neighboring 
W4 contains the most active star-forming sites: W3 North, 
W3 Main, W3 (OH), and AFGL 333. Activity in most of 
these regions might have been triggered by nearby clusters 
and high-mass stars (e.g., lOev et aHl200 5F). More localized 
high-mass star formation is found in the western fields (e.g., 
the KR 140 H II region), which show indications of a more 
quiescent or isolated evolution with spor adic or sequential pe- 
riods of star formation (e.gjRivera-Ingraham et all 201 ll see 


IR i ve ra-1 n g rah ami 20121.1 Rive ra-1 n g rah am et al.12013 (Paper I), 

and lMegeath et al.l2008l for a detailed description of the cloud 
and a review of recent literature). The most prominent re¬ 
gions in W3 have been labeled in Figure Q] which shows the 
column density map at a resolution of ~ 36" produced with 
Herschel^] data (Paper I). 

In Paper I we used the column density map to investigate the 
nature of the most massive and highest column density struc¬ 
tures. We found these structures to be most likely the result of 
feedback by high-mass stars, the combined, convergent effect 
of which might ultimately be responsible for the formation of 
the most massive Trapezium-like systems and clusters. 

In this paper we present the second part o f our Her¬ 
schel a nal ysis o f the W3 GMC (iRivera-Ingrahaml 120121 : 
iRivera-Ingraham et al.ll2013l) . Results below support the con¬ 
clusion from Paper I that local stellar feedback appears to be 


19 Herschel is an ESA space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with important 
participation from NASA and the CSA. 
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Figure 1 . Ah 2 column density map of the W3 GMC corrected for contributions from foreground/background material. Labels mark prominent features in the 
cloud. The white cross marks the intersection of the four subfields in W3; counterclockwise from upper left: W3 Main/(OH), AFGL 333, KR 140, and W3 NW. 
The position of VES 735 is marked with a white X. The HDL is the dense region at the east comprising W3 Main, W3 (OH), and the AFGL 333 Ridge. Contours 
supplement the colorbar at the highest column densities, at Ah, ~ [30,60,200] X 10 21 crrT 2 . 


a major player not only in the formation of clusters of high- 
mass stars but also in the overall process determining the char¬ 
acteristics and local evolution within a GMC. 

In Section [2] we provide a brief introduction to the Her- 
schel datasets and analysis techniques. Section [3] introduces 
the probability density functions (PDFs) and cumulative mass 
distributions (CMDs) and the procedure for selecting the 
structures associated with star formation. Appendix [A] dis¬ 
cusses the effects of foreground/background (non-GMC) ma¬ 
terial on the analysis and Appendix IB1 describes some details 
of the interrelationship between the PDFs and CMDs used. 
Section |4] embarks on the interpretation of structure. The re¬ 
sults of our comparative analysis of the subfields in W3 are 
included in Section [5] Section [6]presents new evidence from 
the PDFs that could be used to trace and constrain the birth¬ 
places of clusters of high-mass stars in a given region. We 
conclude in Section [7] with a summary of the key findings re¬ 
lated to cloud structure and star formation properties of this 
cloud and how they relate to the results presented in our pre¬ 
vious Herschel study of W3 regarding the high-mass star for¬ 
mation process. 

2. DATA PROCESSING AND MAPS OF COLUMN DENSITY 

The W3 GMC was observed with Herschel dPilbratt et al.l 
120101) as part of the HOBYS0 Key Programme (He rschel 
imaging survey of OB Young Stellar objects; iMotte et all 


120101) . In Paper I we presented the Herschel observa¬ 
tions (SPIRE/PACS parallel scan ObsIDs: 1342216019, 
1342216020; bright mode SPIRE ObsIDs (for saturation cor¬ 
rection): 1342239797, 1342239796), ancillary data, and the 
techniques adopted to create and analyze maps of dust optical 
depth and temperature for this cloud. We recall that these 
maps were made by assuming a constant dust temperature 
along the line of sight and fitting spectral energy distributions 
(SEDs) pixel-by-pixel using the Herschel dust emission maps 
at wavelengths >160 fim convolved to the resolution of the 
500 yum map (~36"; ~0.35pc at a distance of 2kpc). To de¬ 
scribe the SED we assumed a fixed dust emissivity index of 
P = 2. 

The optical depth and dust temperature maps were also 
corrected for the effects of emission from dust in the fore¬ 
ground/background of W3. This process removed the non- 
GMC components to reveal more closely the properties of 
the interstellar medium within the GMC. The correction was 
accomplished by estimating the contribution to the emis¬ 
sion I v at each Herschel band from dust traced by atomic 
and molecular gas (H I and CO emission, respectively); see 
iRivera-Ingrahaml d2012h and Appendix B of Paper I for a de¬ 
tailed description of the steps and assumptions associated with 
this technique. This work uses these foreground/background 
interstellar medium (ISM)-corrected maps as the default for 
our analysis. 
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2.1. Alternative Representations of Column Density 

To convert the observable, the dust optical depth, into gas 
column density, A'h 2 (Figure [TJ, we assumed a dust opacity 
of 0.1 cm 2 gm 1 at 1 THz and a mean atomic weight per 
molecule of /./ = 2.33. The latter may be as high as 2.8 
( Roy et aP 120131) but this slight inconsistency was retained 
here so that the column densities could be directly compared 
with results in other Hersch el and HOBYS fie l ds based on 
the same assumptions (e.g.. lAndre et al.ll2010t iMotte et akl 
IMolilHill et alJl2()TTI:lHin et al.ll2012l) . While f3 is not defini¬ 
tively known in molecular clouds and recent work poin ts to 
o pacity variatio ns that should be considered (iRov et al.1120131 
iPlanck Collaboration XIil2014l) . none of these systematic un¬ 
certainties in scale should impact our results significantly, nor 
is such an investigation within the scope of this paper. 

Many results in the literature relevant to the topic of this 
paper are given in terms of magnitudes of dust extinction Ay 
although that is not usually directly observable. For the pur¬ 
poses here we adopted the common approach of transforming 
1 Vh 2 to A v using Ah, = 0.94 x 10 21 Ay cm -2 , although this is 
calibrated only for lower column densities and largely atomic 
lines of sight (iBohlin et al... 197 8). Converting directly from 
submillimeter dust optical depth to Ay would seem more de¬ 
sirable/less contrived and indeed the challenges of doing this 
have been discussed (iPlanck Collaboration Xll 120141) . How¬ 
ever, again our conclusions here should not depend on these 
precise details. 

For the corrected maps applying to the GMC material at a 
common distance of 2 kpc, we can convert the column density 
as parameterized by Ay into a mass per pixel: M p = /Ay. 
For the 9" pixels used and the assumptions described above 
(distance to W3 and mean atomic weight per molecule) we 
find / = 0.13 M e mag~ 1 . 

3. METHODOLOGY: IDENTIFICATION AND CHARACTERIZATION 
OF STAR-FORMING STRUCTURES 

To characterize the density structures in W3 we created 
PDFs of the column density maps and from these mass dis¬ 
tributions. 



1 10 100 
Av 


Figure 2. Probability Density Functions (PDFs) in W3. The logarithmized 
version binned in log Ay is plotted. Top: PDFs for each of four fields in W3. 
Errors bars are shown here as a reference. Bottom: Global PDF for the entire 
W3 GMC with a fitted lognormal function (red) and power-law tail. 

PDFs quantify the fraction (or the probability p) of mate¬ 
rial in the cloud having a column density in the range A/, 


and A/ 2 + AA'n,. PDFs have been used extensively as an¬ 
alytical tools for describing the distribution of mass in re¬ 
gions of low-mass star format ion (e.g., iFroebrich & Rowlesl 
2010; Kai nulainen et al.1 120091) and in regio ns of high-mass 
star formation as well (e.g., HOBYS studies: iHill et al .||20 fTt 

ISchneider et al.ll2012HHill et al.ll2012tlRusseil et al .11201 3l). 

The lower panel of Figure [2] shows the PDF for the entire 
W3 GMC. What is plotted is a “logarithmized” version of the 
PDF, P = Pf io, where the independent variable (for the bin¬ 
ning) is now log Ay (base 10). 

The shape of PDFs is not always straightforward to in¬ 
terpret, as various physical processes imprint on its struc¬ 
ture. For example, interstellar turbulence likely determines 
the lognormal distribution of low_column-densities (see the- 
oretical work of e.g.. iKlessenl l2000l and observations of 
iKainulainen et al.1 120091 for all cloud types: quiescent, low- 
mass, and high-mass star-forming clouds). Following the ap¬ 
proach used in previous studies, the main peak in the above 
PDF and those below was fitted with a lognormal distribution. 
All fits were carried out using a non-linear l east-squ ares mini¬ 
mization idl routine based on mpfit ( I M ark wardtH2009l) and the 
result plotted along with the underlying data. 

External compressio n can cause a broadening of th e PDF, as 
predicted in models of i Fede rrath & Klessen i ( 12013!) and seen 
in Orion B (ISchneide r et al.||2013b and i n clouds associated 
with H II regions dTremblin et al.ll2014l) . The shape of the 
PDF might therefore be used to infer concrete basic physical 
properties of the cloud and individual structures whose mate¬ 
rial is traced by the PDF Ie.g.. iFischerall2014l) . 

Of all processes, gravity plays the most important role at 
higher column densities and for low-mass star-forming re¬ 
gions causes a clearly de f ined power-law like tail in the PDF 
(IKainulainen et all 120091 I Andre et~akll201 It ISchneider et al.l 
1201 3h . High-mass star-forming regions often also show tails 
at high extinction that can also be a ccurat ely modelled as 
power-laws (NGC6334; iRusseil et al.ll2013h . Such high-Ay 
tails can also have, however, more complex shapes as well, 
as has been found here within W3 and in other fields (see 
also lHill et al.l201 ltlSchneider et al.l2012l) . Indeed, processes 
such as turbulence, gravity, feedback/compression, magnetic 
fields, intermittency of density fluctuations, a non-isothermal 
gas phase, properties of the cloud formation processes, and 
even line-of-sight effects, could lead to complex substructure 
in the PDF tails, such as “breaks” and peaks, that might not 
necessarily be well represented with a simple power-law func¬ 
tion. In the next sections, we will discuss in more detail 
how stellar feedback might also imprint on the PDF shape. 
Effects on the PDFs produced by contamination from fore¬ 
ground/background material along the line of sight, avoided 
here, are described briefly in AnncndixfAl 

Column densities theorized in other studies to be associ¬ 
ated with high-mass star formation are very high, for ex- 
ample X = 0.7gem 2 to produce a star with M ~ 10M 0 ; 
iKrumholz & McKeel2008l. which corresponds to Ah. ~ 180x 
10 21 cm -2 or Ay ~ 200 in our maps and PDFs therefrom. 
Cloud column densities of this order, typically observed only 
towards regions of active high-mass star formation, are rare 
and therefore there are often poor statistics at the high extinc¬ 
tion end of PDFs. In our explorations we found that the cumu¬ 
lative form of the mass distribution, the CMD £i_.e., the total 
mass above any given magnitude; e.g., IFroebrich & Rowlesl 
120101) . allows for a complementary and often more straight¬ 
forward analysis of the higher end of the PDFs. 
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Figure 3. Cumulative mass distributions for the W3 field (black: solid line) 
and for each of the four fields: W3 Main/(OH) (black: dotted line), AFGL 
333 (red), KR 140 (green), and W3 NW (blue). The column densities Ah? 
can be compared directly with Figure |T| In terms of Ay, the x-axis range is 
very similar, 0 to 320 mag. 



T [K] 


Figure 4. Like Figure[2]but for the logarithmized PDF of dust temperature. 
Top: for each held in W3. Bottom: global PDF of the W3 GMC. 

The large differences between neighboring regions in W3 
(e.g., the dense and active HD L vs. the more quiescent and 
difFuse western fields; iRivera-Ingraham et al.ll20 lTk Paper I) 
make it necessary to quantify how in-cloud local conditions 
affect the star formation process. To this end, we have carried 
out an analysis of individual areas within W3 itself, these be¬ 
ing the four fields into which we divided this cloud in Paper I: 
the W3 Main/(OH), AFGL 333, KR 140, and W3 NW fields 
(Figure [I]). 

The PDFs of each of the four fields in W3 are given in the 
upper panel of Figure |2j the characteristics of the global PDF 
of W3 in the lower panel clearly depend on the contributions 
from each of the four fields. 

Figure0shows the corresponding CMDs, very different for 
the four fields. In this figure this is expressed in terms of 
Ah, for more direct comparison with the high column density 
regions in Figure Q] 

4. PDF INTERPRETATION AND CLUES TO CLOUD STRUCTURE 

4.1. Analysis and Characterization of PDF Substructure 

While various factors can influence the shape of a PDF, 
the column density above which it deviates from a lognor¬ 


mal distribution into a power-law (i.e., the break) has been 
commonly interpreted as the point separating the turbulent 
medium from the star formation regime. The PDF break can 
therefore be understood as the column density above which 
star-forming str uctures start to dominate over the local en- 
vironment fe.g. . iKainulainen et all [20091 : I Andre et~aj] 1201 It 
iBallesteros-Paredes et al.ll201 ll and references therein). 

This break is therefore an essential parameter for constrain¬ 
ing the star formation properties of individual fields. We 
have found that such breaks can be identified and character- 
ized by fitting linear fu nctions to the log of the CMD. See 
iRivera-IngrahamlfeO 1 2l and Appendix [B] for a more detailed 
description of the concepts and methodology used. Although 
the accuracy of the estimate of the break can be affected by the 
binning (A,4 v = 0.5 mag) and fitting procedure (uncertainties 
in linear slopes are ~ 5 %), a major source of systematic un¬ 
certainty for the lowest breaks will arise from the corrections 
applied to the maps to remove line-of-sight, non-GMC mate¬ 
rial, as explained in AnnendixlAl 

In our application of this approach, the division into sub- 
fields from Paper I was used, as it separates regions with 
dramatically different characteristics. Our choice of subfields 
was determined by unique differences in column density and 
temperature distributions, stellar activity, and stellar content. 
Constraining the origin of these differences is crucial for un¬ 
derstanding the onset of high-mass star activity in the very 
specific regions of the HDL. The temperature differences can 
clearly be observed, for instance, in the temperature PDFs for 
each field in FigurelU The W3 Main/(OH) field is the most ac¬ 
tive and warmest, as well as the only one with ongoing high- 
mass star formation. At the other extreme, the KR 140 and 
W3 NW fields are cold and quiescent. The AFGL 333 has 
intermediate properties, with some high-mass stars, and yet is 
much more quiescent than W3 Main/(OH). 

The derived values and parameters for the linear fits (slope 
and intercept of the lines used to represent the mass distribu¬ 
tions) are given in Table Q] and the fitted lines are shown in 
the right panels of Figures [5] and [6] An expanded view of the 
CMDs in Figure0has been included in Figure^ 

The Figures show that the first breaks (Ay(SF)) derived 
from the mass distributions approach those in the PDFs, while 
the shapes of the components can effectively account for vi¬ 
sual changes (bumps) in the tails of the PDFs. The accuracy 
for determining the break clearly depends on the sharpness 
of the transition between the different (linear) regimes in the 
mass distributions (i.e., those fitted with linear functions). A 
slow transition between two regimes (i.e., not sharp, but pro¬ 
gressive and occurring over a broader extinction range, as ob¬ 
served for those in the HDL; Figure[7]) would result in a higher 
degree of uncertainty, as the identification of the end and start¬ 
ing points of the linear regimes becomes less clear. A higher 
degree of uncertainty for the first break of the KR 140 field is 
expected due to this effect and the fact that this field is also the 
one with the most severe line-of-sight (molecular) contamina¬ 
tion. 

Prominent differences can be seen between the PDFs of 
the two regions in the HDL and those of the western fields. 
W3 NW and KR 140 (see Figure^ show simpler, well-fitted 
PDFs. While a classical power-law function would not be 
able to reproduce the PDF tails of the western fields, our fits 
derived from the linear functions to the mass distributions can 
reproduce them well. The W3 Main/(OH) and AFGL 333 
fields (see Figure[5]) have much more complex tails extending 
to high extinction. While our fits do not manage to repro- 
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Table 1 

Breaks and Parameters 2 for Linear Fits to the CMDs 


Field 

Mass b 
(10 4 M 0 ) 

Breakl c 

(mag) 

Break2 d 

(mag) 

Slope/Intercept (1) 

Slope/Intercept (2) 

Slope/Intercept (3) 

Main/(OH) 

6.0 

13.0 

38.5 

-0.0469/+4.8823 

-0.0154/+4.4656 

-0.0032/+3.9958 

AFGL 333 

6.3 

6.0 

22.5 

-0.0920/+5.0419 

-0.0360/+4.7095 

-0.0084/+4.0903 

KR 140 

5.9 

6.5 


-0.2120/+5.3082 

-0.0945/+4.5196 


W3NW 

4.8 

6.0 


-0.1465/+5.0572 

-0.0844/+4.6769 



2 Derived from log M = a Ay + b, with M in M 0 . 

b Total (foreground/background ISM-corrected) mass of GMC in this field (see also Paper I). 
c Ay(SF) rounded to nearest 0.5. 
d Ay(HTB) rounded to nearest 0.5. 




a W3 Main/(OH) field 




b AFGL 333 field 


Figure 5. PDFs (left) and CMDs (right) for the two fields in the (eastern) FIDL. CMDs: Solid green, red and blue lines are the best linear fits to the data. 
Vertical dash-dotted lines mark the breaks in the distribution (Ay(SF) and Ay(FITB); see text). PDFs: Vertical dash-dotted lines and green, red, and blue curves 
are those derived from the CMD (Appendix 151. Magenta curve is the best-fit lognormal function to the PDF peak. Solid vertical lines are Ay(FlB), where the 
flatter components start to dominate. 


duce this complexity as accurately as those for the western 
fields, they successfully locate the Ay(SF) break and man¬ 
age to trace the overall variation of the tails with increasing 
extinction. Note that, like for the western fields, the tails of 
these PDFs would not be well fitted with a classic power-law 
function either. 

The fits to the mass distributions reveal the presence of a 
second major break in the PDFs of the HDL fields (defined 
in this work as Ay(HTB), for “High-extinction Transition 
Break”)- This break marks the transition point between the 
classical power-law tail of the PDF and a “flatter” regime, fit¬ 
ted in the mass distributions with a third linear function (blue 
slopes; Figure[5]i. As with the first break in the KR 140 field. 


this second transition was found to be considerably wide and 
progressive (i.e., characterized by a smooth, slow transition) 
for both HDL fields, therefore resulting in a higher degree of 
uncertainty when determining the true transition point. The 
regions traced by these intermediate (transition) extinctions 
are identified with the structures and shells in W3 Main/(OH) 
and the AFGL 333 fields, which host the most massive ridges 
and clumps in the W3 GMC. 

The actual break at which the flat regime at high extinctions 
starts to dominate in the PDF (Ay(HB)) is marked by the point 
where the mass distribution is properly described by the third 
linear function, whose parameters have been included in Table 
IH This extinction is also coincident with the point where the 
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A v 

a KR 140 field 


A v 




b W3 NW field 

Figure 6. Same as Figure [3] but for the two western fields. These two fields lack the flattening of the PDF at high extinction that characterizes the HDL fields. 



Figure 7. Expanded view of the low extinction range of the CMDs in Figure [3] highlighting the breaks and the best linear fits to the data for W3 Main/(OH) 
(left), and AFGL 333 (right). 


two linear fits (blue and red) in the mass distributions intersect 
in the PDF (Figure^. These points are Ay(HB)~ 45 mag and 
95 mag for AFGL 333 and W3 Main/(OH), respectively. The 
presence and nature of this possible break is important due 
to its link with high column density material, including high- 
mass star-forming regions. 


Kolmogorov-Smirnov (KS) probability tests). Considering 
also their dramatic differences in (high-mass) stellar activity, 
characterization of the different Ay regions, with the aid of 
the PDFs, is essential for understanding the processes driving 
current star formation and the history that led to these (suc¬ 
cessful) conditions. Below, we discuss each Ay range in turn. 


4.2. Interpretation and Clues to Cloud Structure 

While all fields have comparable total mass, those in 
the HDL have distributions of material at higher extinc¬ 
tions that distinguish them, not only from the KR 140 and 
W3 NW fields, but also between themselves (e.g., based on 


4.2.1. Low Ay Range: Constraining Environmental Conditions 

Compared with other high-m ass sta r-forming regions (e.g.. 
Rosette cloud; ISchneider et al.l 120121) . the properties of the 
PDFs/mass distributions of W3 indicate that material with ex¬ 
tinction Ay ~ 3 mag comprises a typical GMC environment. 
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Table 2 

Characteristics of On-going and Potential 
Star-forming Structures 


M 

(10 3 M 0 ) 

% a 

<T> 

(K) 

Area 

(pc 2 ) 

% 

W3 Main/(OH) 




b 20.5 ± 0.12 

34.2 

21.1 

39.8 

7.6 

c 4.9 ± 0.04 

8.2 

18.3 

11.1 

2.1 

AFGL333 





34.1 ±0.09 

53.9 

17.1 

151.1 

24.7 

19.8 ±0.07 

31.3 

15.5 

80.4 

13.1 

KR140 





9.5 ± 0.03 

16.2 

15.5 

59.4 

6.6 

7.3 ±0.03 

12.4 

15.1 

42.9 

4.7 

W3NW 





15.9 ± 0.04 

32.8 

15.0 

97.6 

15.2 

14.9 ± 0.04 

30.7 

14.8 

89.4 

13.9 


a % of total field. 
b A v > A V (SF) (Table [T}. 
c A v > A V (SF) and T < T sm . 


or common plateau, on which the structures with star-forming 
potential are observed. Note that the critic al extinction for 
core formation of A v ~ 8 mag (1 Andre et al.ll2010t) relies on a 
filamentary environment, so the above estimate represents the 
environment of the filament itself. 

Excluding W3 Main/(OH), the only field with confirmed 
on-going high-mass star formation (e.g., HC H II regions), 
the actual transition point to the gravity-dominated regime in 
W3 is found to be <Ay(SF)>^ 6 mag. M ore than 9 0% o f 
the Class 0/1 YSO population (IRivera-Ingraham et al.ll201 11 
is contained above the Ay(SF) of a given field (except in KR 
140, where ~ 30% of the population is below the Ay(SF)). 
The PDF breaks are therefore suitable limits for identifying 
the major current and potential star-forming sites in the GMC, 
separating these from typical environmental column densities 
in the cloud. The amount of star-forming cloud mass in each 
field, based on these limits, is shown in Table[2] 

Similarly, we can also define typical environmental tem¬ 
peratures, or T env , for each field. The temperatures associated 
with material below the Ay(SF) coincide with the peaks of the 
temperature PDFs shown in Figured] i.e., T 19.5 K, 17.5 K, 
16.0 K, 16.0 K for W3 Main/(OH), AFGL 333, KR 140, and 
W3 NW, respectively. Note, however, that a key property of 
the youngest high-mass star-forming regions in W3 is their 
association with high-column density material with tempera¬ 
tures highly above environmental (Paper I). Exclusion of such 
warm material when identifying young regions would clearly 
underestimate the true star-forming potential of the GMC, es¬ 
pecially in terms of high-mass stars. 

4.2.2. High A v Range: Searching for Sites of High-Mass Star 
Formation 

Characterizing the nature of the second break at high Ay in 
the PDFs is crucial due to its possible association with regions 
of high-mas s star fo rmatio n. While similar breaks were also 
observed bv lSchneider et all (120121) for the Rosette field con¬ 
taining the highest column densities (reaching Ay ~ 70 mag 
in their maps), the AFGL 333 and W3 Main/(OH) fields reach 
extinctions 2-6 times this limit, and second breaks about 2.5- 
4.5 that observed in Rosette. In the following sections, we 
suggest that these second breaks originate from, or are in¬ 


fluenced by, the effects of an external dynamic process. In 
the case of W3, this process is dominated by stellar feed¬ 
back and associated triggering (compression) from high-mass 
stars, acting on already relatively dense star-forming struc¬ 
tures (e.g., shells). Therefore, the structure of the PDFs of W3 
provide additional evidence for the convergence constructive 
feedback scenario of high-mass star formation (Paper I). 

5. DISCUSSION: FAVORABLE CONDITIONS FOR HIGH-MASS STAR 
FORMATION 

The W3 GMC offers a unique opportunity to investigate 
star formation under different environmental conditions, from 
high density structures with high stellar feedback (eastern 
HDL), to more diffuse and quiesc ent (western) regions with 
localized star formation (iRivera-Ingraham et al.ll20lH) . In the 
following sections we compare the fields in W3 to constrain 
the conditions that have led to the (rare) onset of high-mass 
star formation in this GMC. 

5.1. The Environmental Factor Ay( SF) and the Role of 
Stellar Feedback in Determining Local (in-Cloud) and 
Global Evolution 

The HDL, hosting the only young high-mass population in 
W3, comprises up to ~ 70 % of the dense (above environmen¬ 
tal limit; Ay(SF)) material in the GMC. For comparison, W3 
Main/(OH) and AFGL 333 contain up to ~ 2.5 times more 
mass for potential star formation than the neighboring west¬ 
ern fields KR 140 and W3 NW. 

Assuming that a third of the total mass of W 3 wi th 
Ay >A V (SF) is transformed into stars (lAlves et al.l 12007 ). 
the HDL would therefore have a total mass fraction involved 
in star formation (MasssF) or “MSF” in [Froebrich & Rowlesl 
1201 Oh of ~ 15 % at a resolution of ~ 0.35 pc (c.f. the west¬ 
ern fields: Masssv ~ 7.5% of their total mass M tot ~ 1.1 x 
10 5 M Q ). The W3 GMC as a whole shows a MasssF ~ 11.5 %. 

While the mass of the W3 GMC is comparable to that of 
Auriga 1 or Cepheus, these clou ds only have Mass.sF of jus t 
0.19% and 0.26%, respectively dRuebrich_&Row|es[2010b. 
The maximum MasssF found by iFroebrich & Rowies (12010b 
(in their cloud sample) is ~ 10% (Corona Australis). Com¬ 
pared to these low-mass star-forming clouds, W3 appears to 
have an anomalously high proportion of mass involved in star 
formation, with a very high potential to form new stars in the 
next 10 6 yrs (IFroebrich & Rowlesl l2010h despite its already 
significant ongoing star activity. 

W4 shows signatures of prominent stella r activ ity that has 
been ongoing for at least ~ 6 — 20Myr (lOev et aO2005b . Re¬ 
gardless of the distribution of material in the W3 region prior 
to the formation of W4, the scenario of successive episodes of 
(high-mass) star formation and bubble development described 
by these authors suggests that W4 not only influenced the star 
formation process in W3 at the distance of the HDL, but also 
that the HDL, the densest structure in W3, originated due to 
this same activity of bubble/shell expansion, redistribution, 
and compression of material. The idea of a triggered origin for 
the HDL (e.g., by compr essio n) has already been suggested 
in previous studies (e.g.. iMoore et al.ll2007b . and is strongly 
supported by extensive observational evidence. Its location, 
parallel to the W4 H II region, the morphological and phys¬ 
ical characteristics of the HDL presented in this work, as 
well as the W4-W3 s t ellar age and cluster distribut ions (e.g., 
ICarpenter et al.ll2C)OOt IRivera-Ingraham et al.ll20 1 T[) . indicate 
a clear influence of W4 on the material as well as the stellar 
activity in the HDL over an extended period of time. Based 





































8 


Rivera-Ingraham et al. 


Table 3 

Average Environmental Properties of YSOs a 
in the W3 GMC. 


YSOs 

<T d > 

(K) 

< Ay > 

(mag) 

W3 Main/(OH) 

+ AFGL 333 


All types 

18.1 

25.6 

Class 0/1 

16.9 

56.6 

Class II 

18.5 

14.8 

KR 140 + W3 NW 


All types 

15.3 

9.1 

Class 0/1 

14.8 

12.9 

Class II 

15.5 

7.9 

a Catalog 

1, all 

flags from 


IRivera-Ingraham et all 1201111 . 


on this result, and with the HDL dominating the contribution 
of dense material in the GMC, stellar feedback can therefore 
explain the high MasssF in the HDL and therefore in W3 as 
a whole. Feedback as the driver distinguishing the structures 
in the HDL from those in t he western fields is also supported 
by molecular observations (iPoIvchroni et al.ll2012 h. The key 
issue remains as to whether or not feedback can also explain 
the local differences between fields, and why high-mass star 
formation is exclusive to just some particular regions in the 
HDL (W3 Main/(OH)). 

Theoretical models and simulations predict that stellar feed¬ 
back can indeed have a significant impact on the star forma¬ 
tion process. Depending on environmental conditions and the 
location and number of high-mass stars, feedback strongly 
disturbs cloud morphology, shifting and redistributing mate¬ 
rial and creating new dense regions potentiall y suitab le for 
low a nd high-mas s star formation ali ke (e.g .. (Whitworth et al.l 
fl99i IPal^~etaD 120071: IWalch et al.ll20l3h . It can also af¬ 
fect the average stellar mass, by dispersing or concentrat¬ 
ing the local material needed for accretion, while increas- 
ing the total number of young stars of a given region (e.g., 
iFederrath et alJl20l4t IDale et all 120151) . Its effects on local 
star formation from winds and (especially) ionization (e.g., 
IDale et ahll2013l) are, however, dependent on distance and en¬ 
vironmental density, with its destructive effe cts bei ng highly 
minimized in the densest en vironments (e.g. . IDale & Bonnelll 
1201 UlNgoumou et al.ll2015l) . Ultimately, feedback might reg¬ 
ulate the global star formation process at a range of spatial 
scales through the input of turbulence (e.g., iKlessen et all 
12004 Boneberg et al.li20I~5h and a complicated balance be¬ 
tween destructive (e.g., material dispersal) and constructive 
effects (e.g., compression, material accumulation and creation 
of de nse structures; e.g., IDale et all 120071: iKrumholz et al.l 
120141 and references therein). 

Based on our observations and the general theoretical pre¬ 
dictions in the above studies, we suggest that stellar feedback 
from high-mass stars is currently the main factor distinguish¬ 
ing the observed differences in star formation between fields 
in the W3 GMC, disrupting the local environment and altering 
the characteristics, onset, and evolution of the star formation 
process in a given region. This effect can be inferred from 
Tabled which summarises the environmental differences be¬ 
tween W3 YSO populations. Here we used the /V H , and T of 
the pixel in the Herschel maps coincident with the YSO co¬ 
ordinates as representative of the local conditions in which a 
given YSO resides (pixel size~0.1 pc). 

The mean extinction of Class 0/1 and Class II candidates in 


the HDL (IRivera-Ingraham et al.ll20lH) is Aq/i^hdl > 4.4 x 
Ao/i-west and A U - HDL » 1.9 x A U - W e S t, respectively (Table©. 
The triggering process that originated the HDL, in combina¬ 
tion with the activity from the local high-mass stars in W3, 
have therefore provided the denser (and warmer) conditions 
suitable for the onset of the most vigorous and richest stellar 
activity currently observed in the W3 GMC. 

Similarly, comparison of the local differences between 
Class 0/1 and Class II YSOs in a particular region (A< Ay >, 
or the change in local column density with time), can con¬ 
strain possible environmental effects on YSO evolution. As¬ 
suming a total lifetime for the Class 0/1 (+flat SED) and Class 
0/1 + Cl ass II phases of ~ 0.9 Myr and ~ 2.9 Myr, respec¬ 
tively (lEvans et al.ll2009l) . then for co-eval evolution (same 
age) Class II YSOs in the HDL could leave or have their en¬ 
vironment disrupted (e.g., higher external activity) up to ~ 8 
times faster than in the western fields, whose Class 0/1 and 
Class II candidates co-exist in similar (cool) environments and 
comparable column densities. In a more conservative scenario 
in which Class II sources in the HDL are the oldest (2.9 Myr 
old) Class II population in W3 (while those in the western 
fields have just been formed; i.e., 0.9 Myr old), then Class 
II sources in the HDL still dissociate from their primordial 
material ~ 2.5 times faster than those in the western fields. 
According to theoretical models, this apparently negative ef¬ 
fect could result in the new stellar population in the proximity 
of the high-mass stars in the HDL being more numerous, al¬ 
beit with overall lower masses than the more localized stellar 
population formed in the more quiescent western fields. The 
HDL has, however, a significant high-mass star population 
and massive clusters whose origin might also be linked to the 
(in this case, constructive) effects of triggering and external 
events (Paper I). 

Below we summarize the observational evidence from this 
work supporting a feedback-driven model for the local evo¬ 
lution of regions within the W3 GMC. The se re sults aim 
to introduce the basis of a scenario (Section 15.2b in which 
evolution is intimately linked to the balance of construc¬ 
tive/destructive effects of the feedback mechanism, them¬ 
selves highly dependent on local environmental density. 

5.1.1. The HDL: Star Formation in Very High-Density 
Environments 

W3 Main/(OH) is the only field with ongoing (clustered) 
high-mass star formation, despite the fact that the AFGL 333 
is also influenced by the activity in W4. Considering the ef¬ 
fects of stellar feedback on star/structure formation, this in¬ 
tense star formation, as well as the other unique properties 
of the W3 Main/(OH) field, might ultimately be linked to its 
particularly high degree of stellar feedback from current stars 
within W3 on a quite local parsec scale (a few arc minutes on 
Figure 03, including but certainly not exclusively high-mass 
stars in IC 1795 and those powering W3 Main. This appears 
to the dominant mechanism distinguishing the current state 
and evolution of the different fields. The key role of local 
high-mass stars acting already within the W3 complex, rather 
than external activity from W4, for determining the current 
conditions of the star-forming material within dense environ¬ 
ments is based on the following observations: 

1) W3 Main/(OH) has more extreme environmental condi¬ 
tions, with a T env higher than that observed for AFGL 333. 
This would be inconsistent with the activity in W4 being a 
main factor determining the in-cloud state of the HDL, as the 
latter is closer (in projection), and more heavily irradiated by. 
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the high-mass stellar activity in W4. We find that > 95 % of 
the AFGL 333 field is below the 7' cnv of W3 Main/(OH), ex¬ 
cluding the pilla r east of the AFGL 333 Ridge (YSO Group 
7; IRivera-Ingraham et all 1201 71) . and IRAS 02245+6115, an 
H II region associat ed with a B-type star dHughes & Vineri 
119821: IStraizvs & Kazlauskasll2010l) . 

2) W3 Main/(OH) h as a muc h more significant (and on¬ 
going) stellar activity (IRivera-Ingraham et al.ll201 ll) . as well 
as a higher disruption of the YSO environment. About ~ 80 % 
of Class II sources with Ay < Ay(SF) also have T > T em in 
this field. 

3) While both HDL fields reach column densities of the 
same order, W3 Main/(OH) also has a PDF break A V (SF) 
twice that of the AFGL 333 field. This higher break selects 
the shell-like structures around IC 1795, shaped by the cen¬ 
tral cluster and its ( parse c-sca le) feedback from local high- 
mass stars (IRivera-Ingraham et al.ll20in) . In addition to grav¬ 
itational effects, this PDF break could therefore be directly 
influenced by the magnitude of external/active processes. 

If dynamic (e.g., feedback) processes are key for anoma¬ 
lously enhancing the amount of star-forming (dense) material 
in W3, shifting the position of the PDF break, and driving 
the high-mass star/cluster formation itself (as suggested in 
Paper I), then this common link could, in addition, explain 
why: i) W3, a high-mass star-forming cloud, has a higher 
MasssF th an any of the low-ma s s star -forming clouds in the 
sample of iFroebrich & Rowlesl (120101) . and ii) fields classi¬ 
fied as high-mass star-forming regions have generally a higher 
Ay(SF) than lo w-mass sta r-forming regions (as observed in, 
e.g. JSchneider et akll2012l) . 

5.1.2. The Western Fields: Star Formation in More Diffuse 
Environments 

In the KR 140 field, material with Ay >Ay(SF) is exclu¬ 
sively related to filamentary and shell-like structures; i.e., the 
Trilobite, the shell of the KR 140 H II, and the West Loop. 
(Figure |T} Each of these structures has a morphology con¬ 
sistent with external influence: Radiative-Driven Implosion 
(RDI: iRivera-Ingraham et akll201 ll) . the shell around the O- 
type star VES 735, and the border of a cavity-like structure 
observed in the CO map from_the Canadian Galactic Plane 
Survey ICGPS; iTavlor et a1.120031) at v * 49.3 km s 1 ,respec¬ 
tively. The triggered-like origin of these structures is sup¬ 
ported by the (asymm e tric) dis tribu tion of their YSO popu¬ 
lation ( iRivera-Ingraham et al.ll20lTI) . and column density and 
temperature profiles (e.g., asymmetric gradients). A profile 
example for the Trilobite can be seen in Figure [8] 

The comparison between the two western fields resembles 
that between AFGL 333 and W3 Main/(OH), albeit in a much 
smaller scale. A possible cluster in the central-southern parts 
of the western fields (M. Rahman, priv. comm.) might have 
led to a period of enhanced feedback in the KR 140 field, re¬ 
sulting in higher feedback and a riche r ste llar content than 
in W3 NW; e.g., emb edded clusters (ICarpe nter et al Jl2000l) . 
various (late) B stars (IVoroshilov et al.lll985l) . and the only 
(confirmed) O-star (VES 735) outside the HDL (according 
to SIMBAeFI). This affected the temperature distribution ac¬ 
cordingly, with KR 140 having a much smaller proportion of 
“cold” material in any given extinction range than the W3 NW 
field. 

The western fields therefore provide observational evidence 
of a radically different evolutionary path with respect to those 

21 http://simbad.u-strasbg.fr/simbad/ 


in the eastern layer. Based on their properties relative to 
those of the HDL, such evolutionary difference could be eas¬ 
ily linked to differences in the level and type of feedback in 
the history of W3 NW and KR 140. Here, structures have not 
benefitted from the primordial dense conditions provided by 
the external feedback that the first episodes of star formation 
in W4 provided for the eastern side of W3. Similarly, they 
also lack major sources of local internal feedback comparable 
to that of the HDL. Indeed, compared to the HDL, the local¬ 
ized star formation and the much lower level of local feedback 
in the western fields could also explain the low T env of KR 140 
and W3 NW as a whole, as well as the lack of a second break 
in their mass distributions. 




Offset [10* orcsec] Offset [10 2 arcsec] 


Figure 8. Profiles 200" (~ 2 pc) long through the column density peak of 
the Trilobite centered on position RA/Dec: 2 h 21 m 5. s 3, +61° 27' 2971. Both 
temperature (dashed line; right scale) and column density (solid line; left 
scale) are shown for different orientations of the profile: in the West-East 
direction (top-left); NW-SE (top-right); North-South (bottom-left); NE-SW 
(bottom-right). Black vertical line marks the coordinate center (solid line). 


5.2. The Herschel View of the Star Formation History and 
Evolution of the W3 GMC 

Arising from this work is compelling observational evi¬ 
dence of the dramatic effects that stellar feedback can have on 
the evolution of a cloud. Initially, dense structures with con¬ 
ditions particularly favorable for new star formation can be 
created by local (in-cloud) or, in the case of the HDL, exter¬ 
nal feedback effects. These initial (past) conditions determine 
the properties of the first and subsequent episodes of star for¬ 
mation and associated local feedback events. The cumulative 
effects of this sequence of events will ultimately determine 
the general state of the region at the current epoch being ob¬ 
served. 

Quantifying the balance between destructive/constructive 
effects of stellar feedback, influencing (opposing or aiding) 
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known physical processes associated with star formation, 
such as self-gravity, accretion, and turbulence, is therefore 
fundamental for understanding the star formation process it¬ 
self, from cloud to galactic scales. Our Herschel -based results 
can be used as observational constraints to construct a coher¬ 
ent evolutionary model of the W3 GMC, and therefore of the 
conditions and processes that ultimately originate the exotic 
high-mass stars and clusters. 

Past studies indicate that star formation in the western fields 
must have been initiated at the same time a s (if not before) 
the HDL (e.g.. iRivera-Ingraham et all 1201 1). Results from 
the present work suggest, in addition, that the main difference 
between the evolution and properties of the structures and 
stellar population of the different fields might be the magni¬ 
tude of the local stellar feedback occurring within those fields. 
This difference in feedback level can, however, itself be linked 
to the primordial amount of mass and mean density in these 
fields, and therefore the initial level of star-forming activity. 

A lower initial star formation efficiency (SFE) in a rel¬ 
atively low density medium would result in lower efficien¬ 
cies for inducing star formation in secondary events, and a 
lower likelihood in compressing an already low column den¬ 
sity medium to create structures dense enough for high-mass 
star formation. While this local compression might still result 
in few, relatively isolated, high- mass stars (e .g., KR 140) and 
some localized star formation (iRivera-Ingraham et al.ll201 lh . 
the limited secondary star formation and low environmental 
column densities ultimately restricted the star formation ca¬ 
pability of the western fields to a relatively low level. 

The situation in KR 140 and W3 NW differs from the self- 
enhanced process in regions with initially enhanced column 
densities, like the HDL. Such a set-up would promote an ini¬ 
tially high SFE, resulting in more intense (constructive) com¬ 
pression and triggering events. When acting on an already 
dense environment, this process would lead to the formation 
of new high-ma ss star s, as well as a richer stellar population 
(e.g. JPale et al.|[2007l) . In turn, this population would then be 
more effective in further compressing and confining nearby 
material and therefore further enhancing the star formation 
activity. 

Our c onclusions would agree with the results of 
ICarpenter et al . d2000t) . who found that embedded clusters in 
the W3 region are preferentially located in triggered regions. 
The evidence presented in this work also r esembles the “fire¬ 
works hypothesis” presented bv lKoenig et al.l d2012l) . In that 
model, the stellar feedback by high-mass stars triggers the for¬ 
mation of richer stellar populations, a self-propagating mech¬ 
anism that can spread through the formation of new gener¬ 
ations of high-mass stars when dense material is available. 
In W3, we suggest that this process occurs not only due to 
the original prime conditions created by W4 (i.e., the HDL), 
but also due to the generations of high-mass stars themselves, 
which continue creating the conditions needed for the forma¬ 
tion of new populations of low-mass and, especially, high- 
mass stars and massive clusters. 

The mechanism that we invoke for high-mass star formation 
in dense environments, and more particularly, for W3 Main, 
was first introduced in Paper I. In our proposed scenario, stel¬ 
lar feedback creates new dense material, but the properties 
of the high-mass star population and subsequent triggering 
events, taking place within an already triggered and therefore 
density-enhanced region, ensures that the new dense struc¬ 
tures are also (high-mass) star-forming ones. This is achieved 
in part due to the active process of mass assembly and confine¬ 


ment of the triggering mechanism, that allows for the contin¬ 
uation of the star-forming process. We note that this differs 
dramatically from_the scenarios presented in theoretical mod¬ 
els (e.g., iDale et all 120151) . In these models feedback by O 
stars creates dense material, but this material is incapable of 
forming stars efficiently because it is expelled from the poten¬ 
tial wells that facilitate collapse. Similarly, stellar feedback 
has been predicted in some cases to lower the mass of the new 
stellar population by dispersing local material and disturbing 
accretion, (e.g.. iFederrath et akll2014t) . While more detailed 
simulations are required to test the specific scenario described 
here at sub-parsec scales, our observations suggest that the 
triggering conditions and the dense environments can enhance 
the availability of material and aid the accretion process. The 
combined effects would ultimately lead to the unique popula¬ 
tion of high-mass stars and clusters in W3 Main. 

Star formation in W3 Main start ed and prog resse d inde¬ 
pendently from other regions (iFeigelson & Townslevl 120081 : 
IRivera-Ingraham et al.l 1201 lh . In the former, star formation 
was subsequently enhanced by the on-going low-mass activ¬ 
ity and the local (but large scale - several parsec) triggering 
effect from IC 1795 (Figure [D age ~ 3 - 5 Myr; lOev et all 
l2005llRoccatagliata et al.ll201 If) , itself created by the original 
superbubble activity in W4 (e.g. JOev et al.ll2005i ). This pro¬ 
cess led to the first generation of high-mass stars in the shell 
around this cluster, which ultimately led to the onset of the 
convergent constructive feedback in W3 Main (Paper I). 

6. Av(HB), THE SECOND BREAK IN THE CMD/PDF: TRACING THE 
ORIGIN OF THE BIRTHPLACES OF CLUSTERS OF HIGH-MASS 
STARS 

A second break in the mass distributions is a property 
unique to the HDL fields. In this case, that of the W3 
Main/(OH) field (Ay(HB)~ 95 mag) is observed to be at an 
Ay twice that of AFGL 333. 

Column densities above the Ay(HB) breaks are associated 
exclusively with the two high-mass star-forming regions W3 
(OH) and W3 Main, and the AFGL 333 Ridge, the only 
one of the three without confirmed high-mass star forma¬ 
tion. While a lack of high -mass stars in AFGL 333 might be 
due to a younger age (e.g.,]sikai et al. 2 0071 [Polvchroni et al.l 
1201 2l) . the fact that the mass in the AFGL 333 Ridge above its 
Ay(HB) covers an area equivalent to the area above the same 
extinction in the three clumps in W3 Main/(OH) combined, 
suggests that the lower number of high-mass stars and their 
farther distance from the forming structure might have led to 
lower densities and a smaller degree of compression and con¬ 
finement. 

We observe that all the structures traced by the Ay(HB) 
contour coincide with those we identified in Paper I as most 
likely associated with a dynamic input of material. When act¬ 
ing on a region with already enhanced column densities (like 
the HDL or the shell around IC 1795), boundary high-mass 
stars can be particularly efficient with sub-parsec triggering 
(e.g., AFGL 333 and W3 Main). If the second break in the 
PDF traces those structures associated with such dynamical 
processes (feedback-dominated in the case of the W3 GMC, 
in addition to any additional gravitational inflow of material 
from the local neighborhood this amount of mass might ul¬ 
timately induce), then this feature in the PDF could act as a 
signpost for locating structures with enough mass at high ex¬ 
tinction for possible high-mass star formation. Note, however, 
that having the potential for forming high-mass stars might 
not necessarily translate into actual high-mass star formation 
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itself, as this might occur only under very specific conditions 
(as discussed in Paper I). Similarly, the actual Ay at which 
the second break occurs might vary from region to region, de¬ 
pending on factors such as local environment, strength, and 
direction of the triggering event (e.g., AFGL 333 and W3 
Main). 

Our results link the presence of a high-extinction break 
in the PDF with the effects of an external (dynamic) effect. 
Based on our observations, stellar feedback appears to be the 
major dynamic process acting within and on W3, which is 
the reason why feedback has been specifically mentioned and 
referred to in our discussion as the major driver in the evo¬ 
lution of W3. This stellar feedback-based constructive pro¬ 
cess could also be applicable to other regions (e.g. JXu et al.l 
120131) . In more general terms, conclusions from this work in 
terms of high-mass star formation requirements would hold 
when stellar feedback is replaced or aided by other external 
events capable of recreating similar conditions. This would 
be in agreement with th e co nclusion from the HOBYS study 
of Vela C bv lHill et al.l (1201 ll) . who suggested that the flat part 
of their PDFs could be the result of constructive large scale 
flows. Like the convergent constructive feedback mechanism 
introduced in Paper I, convergence of flows could indeed also 
satisfy the requirement of an active input of material that in 
Paper I we suggested could be the key to high-mass star and 
cluster formation. A study of the applicability of this scenario 
to other regions within the Galaxy, and its associated observa¬ 
tional evidence (e.g., PDFs), is currently the focus of ongoing 
work (Rivera-Ingraham et al. 2015; in prep.). 

7. CONCLUSIONS 

The W3 GMC offers a unique opportunity to investigate the 
formation process in a variety of environments. In this second 
study of W3 with Herschel HOBYS data, we have aimed to 
create a coherent picture of the evolution of this GMC, ana¬ 
lyze its large-scale properties and structure, and further con¬ 
strain the high-mass star formation process as first described 
in Paper I. This study has been carried out by means of a com¬ 
parative analysis of the fields in W3 based on the properties 
derived from the Herschel column density and temperature 
maps. 

The W3 Main/(OH) and AFGL 333 fields show a second 
break in their mass distributions. This break appears to be 
related to the presence of external dynamic processes acting 
on the observed structures. Since this influence is the sug¬ 
gested major mechanism for forming the most massive clus¬ 
ters of high-mass stars (Paper I), this break could act as an 
effective signpost for identifying regions suitable for possible 
high-mass star formation. The actual location of this break 
will depend on the local environmental conditions. 

While the first break of the PDF is expected to be influenced 
by various factors (e.g., gravity), we have presented evidence 
that dynamic processes such as external feedback can also be 
responsible for altering the location of this break. If such pro¬ 
cesses are major players in both shifting the break in the PDF 
and high-mass star formation itself (Paper I), then it could ex¬ 
plain why high-mass star-forming regions have a tendency to 
have a higher break than low-mass star-forming regions, as 
observed in previous studies. 

The combined evidence provided by the YSO population 
and the Herschel datasets suggest that differences in the pri¬ 
mordial local conditions are key for determining the evolution 
and current structural and stellar properties of each field. A 
high initial surface density, mass, and column density could 


allow for a higher initial SFE. The combination of a high 
SFE acting in an already high density region (like the HDL), 
combined with the properties of triggering as a star forma¬ 
tion process, could result in a self-enhancing process in which 
subsequent triggering events lead to an increase of the very 
structures suitable for further star formation. This picture is 
supported by the anomalously high proportion of star-forming 
material in W3 as traced by Herschel , compared to other low- 
mass star-forming clouds. The same events could then lead 
to a r icher population, a fireworks hypothesis as suggested by 
iKoenig et alj (120121) . as well as more massive stars. The par¬ 
ticularly enhanced local large scale feedback observed for the 
W3 Main/(OH) field could therefore explain why this is the 
only field with significant high-mass star formation. 

The western fields on the other hand show only moderate 
stellar feedback. This state would be commonly associated 
with the quiet evolution of a cloud lacking the atypical con¬ 
ditions provided by W4 and the HDL, which have greatly 
enhanced the star-forming potential of the eastern regions of 
W3. 

The combined effectiveness of feedback and similar dy¬ 
namic processes (e.g., constructive convergence of flows) in 
1) the creation of column density structures suitable for star 
formation (as shown in this work), and 2) star formation it¬ 
self (as suggested in previous studies), could then support the 
scenario where star formation progresses simultaneously with 
the formation of their parent structures. This process would 
translate into an increasing energy output (luminosity) of a 
star-forming structure (e.g., core/clump) as the structure itself 
is assembled (equivalent to a diagonal evolution in the L/M di¬ 
agram). Such a model matches the scenario introduced with 
the convergent constructive feedback process in Paper I. 
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APPENDIX 

A. EFFECTS OF (NON-GMC) FOREGROUND/BACKGROUND 
MATERIAL ON COLUMN DENSITY AND PDFS 

While the use of the column density maps corrected for 
foreground/background material can be properly described 
with a single lognormal distribution, it was observed that the 
use of ISM-uncorrected maps resulted in a secondary peak 
at low extinctions (Ay ~ 2.5 mag) for the two northernmost 
fields (W3 Main/(OH) and W3 NW, and therefore also on the 
global PDF of the GMC). This effect can be observed in the 
uncorrected PDFs shown in Figured 



1 10 100 
A 


Figure 9. Top: Probability Density Functions (PDF) for each field in W3 
without correcting for the foreground/background emission. Bottom: Global 
PDF for the entire W3 GMC with two fitted lognormal functions (green and 
red), and a power-law tail (blue). The uncorrected PDF is characterized by a 
double-peak profile. 


Visual inspection of the column density maps revealed that 
the material traced by this secondary peak was associated with 
diffuse material external to the W3 GMC, i.e., below aver¬ 
age internal environmental conditions in the GMC. Indeed, as 
shown in the text, removal of foreground/background material 
effectively eliminated most of the component associated with 
this first peak. Only a “remnant” of a peak is still observed for 
the W3 Main/(OH) field in Figure [9] This leftover feature is 
expected as our correction did not remove contributions from 
the interstellar medium (ISM) local to W3 (i.e., in the velocity 
range of W3 itself), adding to the shift of pixels towards lower 
extinction values when removing the ISM contribution. This 
example emphasizes the need for a careful selection of the 
area chosen for PDF analysis. Note that this type of “double- 
peaked” PDF does not correspond to the ones observed in the 
vic inity of H II regions (lSchneider et al.ll2012l : lTremblin et all 
12014) where the expansion of ionized gas into the molecu¬ 
lar cloud leads to a compressed layer of gas that shows up in 
the column density PDF as a second peak at higher column 
densities. 

While the effects on the PDF depend on the amount of 
ISM correction required for a particular field, this correc¬ 
tion is most important and dominant at relatively low extinc¬ 
tions. This higher degree of uncertainty should be taken into 
consideration when analyzing those regions/structures with 
Ay < 10 mag. Correcting for line-of-sight material broad¬ 
ens and shifts the main peak of the PDFs to lower extinctions, 


which change from Ay ~ 5 mag to Ay ~ 3 mag. The first 
break of the PDFs is observed to shift by Ay ~ 1 mag, while 
changes to the histograms are essentially negligible (within 
binning accuracy of 0.5 mag) for Ay > 30 mag. 

The amount of correction needed will ultimately vary from 
cloud to cloud, and even region to region. The KR 140 field, 
for instance, suffers from the greatest uncertainties due to 
it being severely affected by considerable foreground mate¬ 
rial traced by CO, therefore requiring the largest correction 
of all fields in W3. Indeed, this field has the largest differ¬ 
ence in total mass before and after ISM correction (a factor 
of 1.6 more mass in the uncorrected maps), and shows the 
largest discrepancies in terms of mass above Ay(SF) (~ 3.5 
less mass in the corrected mass for the same extinction level; 
TableO. The mean difference in mass above Ay(SF) between 
corrected and uncorrected maps for the other fields is ~ 20 % 
of the total corrected mass in each field, with the uncorrected 
maps always having more material than the corrected ones for 
any given extinction limit. This situation provides an upper 
limit to the mass uncertainties, as some correction for line-of- 
sight-material is required. Moreover, when dealing with un¬ 
corrected images, the Ay(SF) for each field should be higher 
than those derived from the ISM-corrected PDFs. Therefore, 
the difference between the total mass above Ay(SF) for the 
uncorrected and corrected maps should be smaller than those 
quoted here. None of these uncertainties affect the conclu¬ 
sions from this study. 

B. CONVERTING BETWEEN CMDS AND PDFS 

Similar to the interpretation of t he b reak from the lognormal 
distribution in the PDF (Section 14.11) . it has been suggested 
that the point at which a break occurs in the CMD at low ex¬ 
tinction separates the turbulen t environment from the gra vity 
dominated structures (Av(SF): lFroebrich & Rowlesll2010t) . 

Linear regimes in the CMDs were first selected as those 
regions with slow varying gradient change, relative to those 
with significantly rapid change that should characterize a po¬ 
tential transition or break. These linear regimes were fitted 
in the CMDs with a line function and ay 2 minimization rou¬ 
tine. The final slopes and intercepts in Table Q] are those of the 
linear fits that best represent the data neighbouring the regions 
of maximum gradient change (closest to the breaks) as well as 
the overall shape of the CMD. By fitting the separate extinc¬ 
tion regimes of the log of the CMDs with straight lines (e.g.. 
Figures [5] and |6j, these breaks can be identified as the points 
where linear fits of adjacent regions intersect. The sharper the 
transition between two regimes (with only a small region of 
curvature joining adjacent linear-like regimes) the better the 
constraint on the value of the break. 

Considering the complexity of the PDF tails and the diffi¬ 
culty of fitting them assuming the typical power-law function, 
in this work we explored using the CMDs and the straight¬ 
forward linear fitting method described above as the primary 
approach for locating the breaks. Because a CMD is just a 
cumulative form derived from a PDF, however, any physical 
break should be recognizable when using either of the two 
methods. 

To check for consistency of breaks derived from the CMDs 
with features observed in the PDFs, we transferred our best 
linear fits of the CMDs (T = log M vs. Ay) to the form of the 
PDFs (linear binning). This transfer was done numerically. 
Quite generally, the fraction of material (expressed here as a 
probability p) in a given region with extinction between Ay i 
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and Ay 2 (average Ay) is 


p(A w ) = 


10 r = - 10 y ' 
f n p Ay Bin ’ 


(Bl) 


where Y, is the plotted log of the CMD at Ay ; , / is the con¬ 
stant relating mass per pixel to extinction (Section^, n p is the 
number of valid pixels in the field, and “Bin” is the linear bin 
size in Ay used in the PDF. This can then be transformed to 
the logarithmized PDFs as displayed above. For each line fit¬ 
ted to the CMDs the corresponding colored curves are plotted 
in the left panels of Figures [5] and [6] 

For the specific case in Table [Hand Figures 0 and [6] where 
log M - a Ay + log M(0) or 


In M = .? Ay + In M(0), (B2) 


using c - In 10 and s - -ca, the underlying PDF generating 
function is 

M(0) i 

P(Ay) = ~t~t — exp(—sAy). (B3) 

Ay fll p 

When cast in the format needed for the logarithmized depen¬ 
dent and independent variables used in Figures[5]oi'[H left, the 
curve has the form log Pc io(log Ay) = [ log M(0)+log(-c 2 a)- 
log(/n p )] +al0 logAv . 
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